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We study the dispersion of the tt resonance in the superconducting state within the projected SO (5) model Away 
from the the commensurate momentum, the propagation of the tt resonance creates phase slips in the superconducting 
order parameter. This frustration effect leads to a strong dressing of the tt resonance and a downwards dispersion 
away from the commensurate wave vector. Based on these results, we argue that the the commensurate resonance 
and incommensurate magnetic fluctuations in the cuprates are continuously connected. 



PACS numbers: 74.72. Bk, 61.12.Bt, 61. 12. Ex 

The emergence of a commensurate neutron resonance 
peak below T c is one of the most striking properties of 
the cuprates. There are many different theoretical ap- 
proaches to the problem. Within the SO (5) theory |Q, 
this resonance peak is interpreted as a collective mode 
describing a rotation from the superconducting (SC) 
state to the commensurate antiferromagnetic (AF) state. 
Since the rotation operator is a triplet particle-particle 
operator, which can only couple to neutrons in a SC state, 
this theory predicts that the intensity of the tt resonance 
mode is proportional to the square of the SC order pa- 
rameter H . This fact naturally explains the temperature 
dependence of the resonance intensity observed in exper- 
iments j|, and anticipated [|| the recently observed de- 
pendence of the resonance intensity on the magnetic field 
A c-axis magnetic field creates non-SC vortex cores, 
therefore it reduces the SC order parameter on the av- 
erage. On the other hand, a magnetic field aligned in 
the ab plane has little effects in reducing the SC order 
parameter. From the dependence of the resonance inten- 
sity on the SC order parameter within the SO (5) theory, 
one would therefore naturally expect a reduction of the 
resonance intensity in the presence of a c axis magnetic 
field. 

More recently, incommensurate magnetic fluctuations 
have also been observed in the YBCO superconductors 
Unlike their counterparts in the LSCO supercon- 
ductors, these incommensurate magnetic fluctuations are 
energy dependent, and seem to merge continuously into 
the tt resonance mode. In fact, this behavior is predicted 
by Furusaki and one of us ||, who applied the SO (5) 
theory to the ladder systems. Using a asymptotically ex- 
act renormalization group method, they found that the tt 
resonance mode in the ladder system has a downward dis- 
persion away from the commensurate wave vector (note 
the remarkable resemblance of Fig. 5 and 6 of reference 
H and the experimental Fig. 3 A of Ref. ||). In fact, in a 
earlier work, Poilblanc, Scalapino and Hanke |l(| found 
similar behavior in a exact diagonalization study of the 
ladder systems. The purpose of this work is to identify 
the physical origin of this downwards dispersion relation 



and construct a unified theory of the commensurate and 
incommensurate magnetic fluctuations in the cuprates. 

To understand the basic physics of the competition be- 
tween the commensurate and incommensurate magnetic 
fluctuations, let us first consider a very simple system of 
doped t — J|| — J± ladder. In the J± » J\\ limit, the 
SC state of the doped ladder system can be viewed as a 
coherent linear superposition of a hole pair and a singlet 
pair on each rung. In this picture, the tt resonance is sim- 
ply a linear superposition of a spin triplet excitation on 
each rung. Because the SC state is a broken symmetry 
state, such a excitation can be either created by a spin 
triplet operator or a triplet particle-particle operator, the 
so called tt operator. Simple second order perturbation 
theory shows that the exchange parameter between the 
triplet and a hole pair is — , which has a negative sign, 
while the exchange parameter between the triplet and the 
singlet pair is ^ , which has a positive sign. The positive 
exchange sign prefers a minimal energy wave vector at 
Qmin = tt, while the negative sign prefers a minimal en- 
ergy wave vector at q m in = 0. Since the SC ground state 
is a linear superposition of both the hole pair and the 
singlet pair, this effect leads to a competition between 
these two band minima. In one dimension, this compe- 
tition always leads to a incommensurate band minimum. 
This picture is a crude caricature of the basic physics 
leading to the incommensurate band minimum for the 
tt resonance. For the full calculations in one dimensions, 
readers are referred to the original paper by Furusaki and 
one of us ||. 

In this paper, we shall generalize this physics to two di- 
mensions and study the propagation of the tt resonance in 
the SC state, within the projected 50(5) model (pSO(5) 
model). (In fact, within the pure SC state, there is no 
essential difference between the SO(5) model and the 
projected 5*0(5) model Q.) Just like the physical pic- 
ture outlined above, the propagation of the tt resonance 
leads to a sign reversal of the SC order parameter be- 
hind it. Therefore, the problem is similar to the problem 
of the propagation of a single hole in the antiferromag- 
netic background, which has been studied extensively in 
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the literature |ll|-|14|. In fact, using basically the same 
approximations, we shall show that the quantum fluc- 
tuations erase the string of sign reversals, but leads to 
a non-trivial correction to the dispersion of the it reso- 
nance. Depending on parameters of the model, the cor- 
rection can give rise to a downward dispersion of the 
7r mode, reaching a minimum at incommensurate wave 
vector. Current experiments only show a downwards dis- 
persion of the 7r resonance. The main prediction of our 
theory is that there will be a upward turn in the disper- 
sion after the minimum is reached. 

We begin with the projected 50(5) model defined on 
a lattice, (using the notations of Ref. 0].) 

H = A s ^2 4 fa)*** ( x ) + Ac 2J n i ( x ) n i ( x ) 

X X 

- J s ^2 n a (x)n a (x') - J c 2J n i (x)n i (x') 

<xx'> <xx'> 

+ V Lab{x)L ab {x r ) (1) 



<xx'~> 



where i = 1, 5, a = 2, 3, 4 and A c = A c — fi where fj, is 
the chemical potential. In 50(5) superspin notation, the 
superspin is defined as 
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where t/j and t a are hard-core boson annihilation opera- 
tors for the hole pair and the magnon. Here one lattice 
site of this effective model corresponds to a plaquette in 
the original C u 02 plane, and we have made a shift of 
the momentum vector for t a bosons by (tt,tt). pSO(5) 
model |jj is constructed by projecting out doubly occu- 
pied configuration from the local 50(5) multiplets. Even 
though some members of the 50(5) multiplets are pro- 
jected out, at A s = A c and J c = 2J S , the mean field 
ground state manifold is still 50(5) symmetric, and AF 
can be smoothly rotated into SC with no energy cost. In 
Ref. jlj, we have shown this model gives a realistic de- 
scription of the global phase diagram of the cuprates and 
many of their physical properties. In that paper, the V- 
term in above Hamiltonian was simply ignored because it 
does not impact the mean field phase diagram if only pure 
AF phase and SC phase are concerned. However, one im- 
portant observation in this paper is that this term plays 
important role on the dispersion of n mode. If we limit 
our discussion on one single magnon in Hilbert space, the 
V-term simply describes the hopping exchange between a 
hole-pair and a magnon. The hopping between a magnon 
and a singlet is described in the term —J s n a (x)n a (x'). 
Since we only discuss the motion of an single magnon, 
all of other matrix elements for these two terms vanish 
except 



VL ab (x)L ab (x')(4(x)t+(x')\0 >) = Vt+(s>)t+(x)\0 > 
-J s n a (x)n a (x')(b+(x)t+(x')\0 >) = - J s b+ (x')t+(x)\0 >, (4) 

where b + (x) creates an singlet at site x. Since both of 
J s and V is positive in this model, we see that these two 
hopping processes have opposite signs. Therefore, in a 
pure SC state, which is a coherent state of singlet and 
hole pair on each site, the hopping of magnon flips the 
phase of the local SC order. More precisely, let's take a 
mean field pure SC state, 

<f> = Y[4>+(x), <j)±{x) = {cos6b + (x)±sin6t+{x))\0> 1 (5) 

X 

which has SC order, < n\{x) >= ^sin(26), on each site. 
Taking the hopping term of a magnon between sites x, x' 
we obtain for V = J s , 

(L a b(x)L ab (x') - n a {x)n a {x'))t+{x')(j) + {x) = t+(x)(f>^ (x'), (6) 

in which the SC order at site x', nx(x'), takes value 
< n\{x) >= — ^sin(26). A schematic diagram is shown 
in Fig.|l] to reflect the hopping of a magnon in the SC 
state. The motion of a magnon creates the phase mis- 
match similar to the spin mismatch caused by the motion 
of a single hole in a AF state. 



+ + + + + + + + 

FIG. 1. Schematic diagram of the hopping of a magnon in 

SC state. The arrow represents a magnon while the the signs 

represent the phase of SC order. 

In the Ising limit, the motion of magnon creates a 
phase mismatched string and feels a linear potential 
which results in an infinite effective magnon mass. How- 
ever, similar to a single hole problem, quantum phase 
fluctuations in the SC state can erase the string effect 
and modify the dispersion relation. We employ the same 
method in ref. Jl3|,[l4| to solve the problem. 

Taking the variational wavefunction in Eq.^, a simple 
mean field calculation gives the hole-pair density 



sin 2 {6) 



J c z 



2J c z 



(J) 



where z is the coordination number on the lattice. We 
introduce a Lagrangian multiplier field A to enforce the 
hard-core boson constraint on average. To obtain the col- 
lective phase mode, we define two new boson operators, 



ot\{x) = sin9b(x) — cosOt^x), 
oi^ix) = cos9b(x) + sin9th{x), 



(8) 
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which satisfy ai(x)\$ >= 0, < a 2 (x) >= 1. The 
quadratic Hamiltonian describing the phase fluctuation 
above the mean field state is therefore given by 



H' = H- < H > 



= -j c z £[(i - 7 (g) 



sin 2 (29) 



70?) W (<?)<*! 0?) 



sin 2 (29) 



j(q)(a 1 (-q)ax(q) + h.c.)] 



(9) 



where 7(g) 



cos(q x )+cos(q y ) 



In the following part of 



paper, we normalize all of the energy scale by taking 
( J s + V)z/2 = 1 and define two additional dimensionless 
parameters, 



V - J s 2J C 

X = ~ 77, y 



J s + v' a (J s + v) 

By standard diagonalization, 

ff' = ^e(^+(#(g), 

Q 

where 

/3(g) = u(q)ai{q) - v(q)a+(-q) 



(10) 



(11) 



(12) 

e(q) = yy/(l - ~f(q)) 2 + sin 2 (29)j( q )(l - j(q)) (13) 

(14) 




1 i- l{q) + ^m l(q) 

v(q) = -sign(^(q))\l--+y . (15) 

The hopping of a single magnon can be described by 
the following Hamiltonian, 

H m = — J s ^2(t+(x)t a (x')b + (x')b(x) + h.c.) 

xx' 

+ V^2(t+(x)t a {x')t+{x')t h (x) + h.c). (16) 

xx' 

Taking the mean field expectations and using relations 
of Eqg we get 

H m = Y,(Eb(Q,0)t+(q)t a (q) + h.c.) 
1 

+ £/(*>«)£(*)*«*(* -<?)> (17) 

kq 



where 



E b (q,9) = (x-cos(29)) 1 (q). (18) 



and 



f(k, q) = S zn(2^)[( 7 (fc - q)u(q) + j(k)v(q))f3(q) 

+ (j(k)u(q) + 7 (fc - q)v(q))p + (-q)}. (19) 



Given the Hamiltonian composed of Eq.|TT| and EqJT^, we 
can calculate the dynamic spin correlation function, 

G[x,t) = -i < $|T[i Q (a;,i)& + (x,t)i+(0,0)6(0,0)]|$ > . (20) 

Taking the mean field value on b and using the self- 
consistent perturbation which sums only noncrossing di- 
agrams, we obtain the following Dyson's equation, 



G(k,Lu) = 



where 



cos 2 (9) 



Eo(k, 0) ~ E q F(k, q, 9)G(k -q,u- e(q)) 



(21) 



F(k, q, 9) = sin 2 {29)\ 1 {k - q)u{q) + 7 (fc)«(g)| 2 . (22) 

In this approximation, the vertex is neglected, which is 
small according to the calculation in Ref. [Q for the 
single hole problem. In our model, the problem is just 
more complicated in the sense that the coupling strength 
and the dispersion of phase mode are quite different and 
depend on the density of hole pairs. We numerically 
solve the above integral equations to obtain the spec- 
trum, A(k,oj) — —^G(k,oj) , and in particular the min- 
imum position of the dispersion in the broad arrange of 
parameters. We first observe is that the bare dispersion 



of magnon, Ef,(k), can be removed in Eq,21 by shifting 
u>. Moreover, if y is large, it is better to rescale e(q) in 
numerical calculation, which can narrow the energy ar- 
range required to find a solution. Combine these two 
observations, the numerical convergence is rather fast. 

Throughout the numerical calculation, a sharp coher- 
ent peak is found in all momentum space. The spectral 
function can be generally written as 

A{k,u,9,y) = cos 2 {9)Z{k,9,y)5{uj - Q(k, 6,y)) + A', (23) 

where A' is the incoherent part and £l(k,9,y) defines the 
energy dispersion. From our numerical result, we find 
that the energy dispersion, J7, can be written as 

Sl(x, k, 9, y) = n (x, 0, 9, y) + {E b (k, 9) - E b (0, 9)) 

-6Sl(k,9,y), (24) 

where the first term is the energy at k = (0, 0), the second 
one is the bare dispersion, and the third is the relative 
energy shift contributed by fluctuation correction. Based 
on the symmetry of our model, 5Q(k, p,y) = 6Q(k, 1 — 
p 7 y)). We find that 5fl(k,9) fits very well a product of 
two separated functions, h(9,y) and g(k), where h(9,y) 
is independent of k. Considering the symmetry of the 
lattice, we find the function g(k) can be well fitted to the 
following form, 



g(k x , k y ) = 0.128(sin 2 (k x ) + sin 2 (k y )) 
+ 0.054:(cos(k x ) - cos(ky)) 2 . 



(25) 
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Thus, the minimum of can be analytically determined 
by plugging Eq.|2^ into Eq.^4[ The global minimum al- 
ways happens along the diagonal (it, tt) direction. How- 
ever, along the (tt, 0) direction, there is a local minimum. 
The global and local minimum positions start to shift 
from (0, 0) at a common critical density pi which is ap- 
proximately determined by the equation 



cos(26i) - O.512/i(0/, y) - x = 0. 



(26) 



In Fig.^J, we show the global minimum position, (k m , fc m ), 
and local minimum, (fc m ,0), as the function of the den- 
sity. In Fig.||, we also show the full dispersion of tt mode 
at the density p — 0.4. 

There is another important density in the problem. 
Within the pSO(5) model, uniform SC state can only 
exist for densities exceeding p c . Ignoring small quantum 
corrections, p c is given by 



r - 5 S + 0.5(1 - x) 
0.5y + r + l 



(27) 



where r — J S /(J S + V) and 5 S = j f+ Vz - Therefore, de- 
pended on the parameters, pi could be larger or smaller 
than p c . For example, at x — 0, y = 0.5 and r = 0.25, 
p c > pi for 5 S < 0.3 and p c < pi for 5 S > 0.3. Therefore, 
the relative sizes between pi and p c could tell us the na- 
ture of the states as doping level is reduced. For p < p c , 
a mix state between AF and SC is obtained. If pi > p Cl 
the lowest energy magnetic fluctuations are commensu- 
rate at the transition, and a uniform mix state between 
AF and SC is obtained. On the other hand, if pi < p Cl 
the lowest magnetic excitations are incommensurate at 
the transition, and a stripe state with alternating AF 
and SC order appears. 

Finally we would like to compare our results with 
two classes of theory of the tt resonance, one based on 
particle-particle picture (HJ|] and the other based on the 
particle hole picture Since all calculations are car- 
ried out in the SC state, where these two channels mix, 
distinctions can only be made meaningfully by comparing 
the final physical predictions. In this work, the inten- 
sity of the peak is determined is cos 2 (9) Z (k , 9) , where 
Z(k, 9) is basically independent of the density parame- 
ter 9. Therefore, the intensity of the magnetic collective 
mode is 



the SC order, and the mode energy is always less than 
the SC pairing gap. The particle-hole based pictures |ll| 
can also explain the downward dispersion, however the 
mode terminates at certain wave vector due to Landau 
damping. On the other hand, our current work predicts a 
minimum in the dispersion. The energy at the minimum 
is the spin gap in the system. 
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COS 2 (9) = 1 - p = 



<ni> 



(28) 



where < n\ > is the SC order parameter. This result 
agrees exactly with the prediction of the particle-particle 
picture |qjlql. The energy of the collective mode is in- 
dependent of the SC energy gap, and is expected to be 
temperature independent. Both these properties are in 
contrast with the predictions of based on the particle- 
hole picture, where the mode intensity is insensitive to 
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FIG. 2. The minimum momentum position k m as the 
function of density p at parameters: x = 0, y = 0.5. Solid 
line reflects (tt, it) direction and dashed line represents (tv, 0) 
direction, (we have made a shift back (0, 0) to (it, it) in all 
figures in this paper) 




kx 

FIG. 3. The full dispersion of 7r mode, E n , as the function 
of k x ,k y at density p = 0.4 with the parameters: x — 0, 
?/ = 0.5, z = 0.5, <5 S = 1. 
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